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ABSTRACT 

We analysed the X-ray spectra of six observations, simultaneously taken with XMM- 
Newton and Rossi X-ray Timing Explorer (RXTE), of the neutron star low- mass 
X-ray binary 4U 1636-53. The observations cover several states of the source, and 
therefore a large range of inferred mass accretion rate. These six observations show a 
broad emission line in the spectrum at around 6.5 keV, likely due to iron. We fitted 
this line with a set of phenomcnological models of a relativistically broadened line, 
plus a model that accounts for relativistically smeared and ionised reflection from the 
accretion disc. The latter model includes the incident emission from both the neutron- 
star surface or boundary layer and the corona that is responsible for the high-energy 
emission in these systems. From the fits with the reflection model we found that 
in four out of the six observations the main contribution to the reflected spectrum 
comes from the neutron-star surface or boundary layer, whereas in the other two 
observations the main contribution to the reflected spectrum comes from the corona. 
We found that the relative contribution of these two components is not correlated to 
the state of the source. From the phenomenological models we found that the iron 
line profile is better described by a symmetric, albeit broad, profile. The width of the 
line cannot be explained only by Compton broadening, and we therefore explored the 
case of rclativistic broadening. We further found that the direct emission from the 
disc, boundary layer, and corona generally evolved in a manner consistent with the 
standard accretion disc model, with the disc and boundary layer becoming hotter and 
the disc moving inwards as the source changed from the hard in to the soft state. The 
iron line, however, did not appear to follow the same trend. 

Key words: X-rays: binaries; stars: neutron; accretion, accretion discs; X-rays: indi- 
vidual: 4U 1636-53 



1 INTRODUCTION 

The X-ray spectrum of low-mass X-ray binaries (LMXBs) 
shows evidence of an accretion disc ext ending down clos e 
to the central object (see e.g. review bv lDone et"ai]|2007t ). 
The disc is assumed to be optically thick and geometrically 
thin, and is characterised by a t hermal spectrum with typ - 
ical temperatures of 0.2-1 keV (jShakura fc Sunvaevlll97of ). 
The thermal disc photons are Compton up-scattered in a 
hot coronal gas surrounding the disc and compact object, 
giving rise to the hard comp onent usually seen in the spec- 
trum (e.g. iBarret et af1l2000] , and references therein). How- 
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ever, a fraction of these Comptonized photons is backseat - 
tered, and they irradiate the accretion disc. After being ab- 
sorbed, reprocessed and re-emitted, these photons produce 
a reflection spectrum consisting of a continuum plus a com- 
plex line spectrum (see review bv lFabian fc RosslhoiOl ). Due 
to the high fluorescence yield and abundance, an impor- 
tant feature of the reflection spectrum is the Fe Ka line 
at ~6-7 keV. As a result of disc rotation and relativistic 
effects near the compact object, the reflection spectrum is 
blurred, with the F e line being asymmetrically broadened 
(|Fabian et al.lll989h . As the disc moves close to the cen- 
tral compact object, the relativistic effects become stronger, 
and consequently, the red wing of the Fe line - mainly set 
by the gravitational redshift at the inner edge of the disc 
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- extends to lower energies. Hence, modelling of the iron 
line profile allows us to measure the size of the inner disc 
radius, which determines other spectral and timing prop- 
erties of the source (e.g., disc temperature, characteristic 
frequency of the variability). If the compact object is a 
black hole (BH), modelling of the Fe line profile in the disc- 
dominated s pectral state provides an estimate of the BH 
spin (see e.g. lMilledl2007l : iRevnolds fe FabiarJl200Sf ). assum- 
ing that the disc is truncated at the innermost stable circular 
orbit (ISCO). In case of a neutron star (NS), a measurement 
of the inner edge of the disc sets directly an upper limit 
on the NS radius, as the inner disc radius has to be larger 
than, or equal to, the equatorial radius of the star (see e.g. 
iBhattacharvva fe Strohmaver|[2007l ; ICackett et alj|2008h . 

Since the first detection of a broad and asymmet- 
ric Fe emission line in the Seyfert-1 galaxy MCG-6-30-15 
(|Tanaka et al.lll995h . many similar lines have been observed 
in BH systems, both in Active G alactic Nuclei (AGN) and 
in X- ray binaries (see reviews by iFabian et al . 2000; Miller 
l2007h . In all cases the line appears to be produced very close 
to the BH. Broad Fe lines are also detected in NS systems 
although they are weaker than in BH systems. With instru- 
ments like Chandra, XMM-Newton and Suzaku, relativis- 
tically broadened Fe em ission lines have be e n observed in 
many NS LMXBs (e.g., Piraino et aD |2007t Cackett et al 
200S ; |Pandel et alj|200a|di Salvo et al.ll2009l ; Icackett et al. 
201C ; see also lNg et aljboid for an almost complete list of 



NS LMXBs with Fe lines detected with XMM-Newton). 

Relativistically broadened Fe lines in NS spectra are 
usually fit using a phenomenological model for either a 
Schwarzschild or (maximum) Kerr metric, where the under- 
lying continuum is mo delled with a variety of direct emission 
components (se e e.g. Bhattacharvva fc Strohmavcr 20071 : 
iLin et all 120071 : IPandel et aD 120081 : ICackett et all l200Sft . 



However, as the Fe emission line is a reflection signature, 
it should be accompanied with a reflection continuum ex- 
tending over a wide bandpass; at the same time, the line 
and its reflection continuum should self-consistently deter- 
mine the ionisa tion balance of the disc and include Compton 
scattering (e.g. [George fe FabiarJll99ll : iRoss fe Fabian|["l993l : 
iRoss et al.lll999F An inappropriate modelling of the under- 
lying continuum, as well as not including any of these pro- 
cesses in the analysis, could lead to incorrect results from the 
line modelling. Therefore, a self-consistent ionised reflection 
model should be used to study the continuum and Fe emis- 
sion line in the spectra of LMXBs. (Note that, depending 
on the quality of the data, reflection and phenomenological 
models could describe the data equally well, even leading to 
reasonable results f or the phenomenological model; see e.g. 
ICackett etai]|2010h Furthermore, the shape of the reflec- 
tion spectrum not only depends on the ionisation state of 
the surface layers of the disc, but also on the spectral shape 
of the emission illuminating the accretion disc. Comptonized 
emission from the (disc) corona is the commonly assumed 
source of irradiation, but in the case of an accreting NS, 
the emission from the NS surface /boundary layer, can be as 
well irradiating the accretion disc. lCackett et all (|2010l ) have 
studied a number of broad Fe emission lines in NS LMXBs 
where they used a blurred reflection model in which the 
blackbody component - used to mimic the emission from the 
boundary layer - wa s assumed to be the incident emission. 
ICackett et all (|2010l ) found that this model fits the data well 



in most cases, supporting the idea that t he boundary laye r 
is indeed illuminating the disc (see also iDAi et alj|20ich . 
In our analysis (below), we test whether a disc corona or 
boundary layer was the most likely source of irradiation. 
We note that relativistic (plus Compton) broadening is not 
the only physical explanation for the width and the profile of 
iron li nes in compact objects. For instance. iTitarchuk et all 
(|2003h proposed that asymmetric line profiles could origi- 
nate from an optically thick flow launched from the disc, 
which expands or contracts at relativistic speeds. We will 
not explore this pos sibility in t he pa per, but refer the reader 
to the discussion in lNg et al.l (|2010h . 

We study the Fe emission line and the reflection spec- 
trum in six observations of the NS-LMXB 4U 1636-53, 
taken with the XMM-Newton and the RXTE satellites si- 
multaneously. Besides the direct emission components, we 
fit these data using both phenomenological models and a 
relativistically-smeared ionised reflection model in which we 
investigated different sources of irradiation. In what follows, 
we first introduce the source 4U 1636-53, focusing on quan- 
tities relevant for this paper. In Section [2] we give details 
on the data used, and the procedure we followed to reduce 
them. In Section [3] we describe the models, and explain the 
specific aspects of the spectral analysis. Here we also sum- 
marise the most important findings. We discuss our results 
and their implication on the standard accretion disc model 
in Section [3] There we also address the impact of the NS 
surface/boundary layer, on the reflection spectrum for the 
different spectral states. Furthermore, in that section we also 
discuss the results from the phenomenological line models. 
In Section [5] we summarise our results. 

1.1 4U 1636-53 



ery m 1974 jWillmore et all 1 19741 ; 
iGiacconi et al.l Il974h . 4U 1636-53 has been observed 
over a wide range of wavelengths. Photometry of the optical 
counterpart (V801 Ara) revealed a short orbital period 
of ~3 .8 hr and a companion s tar with a mass of ~0.4 
M© (|van Paradiis et al.l Il990l ; iGiles et al.1 120021 ). X-ray 
studies of 4U 1636-53 have shown a variety of rapid time 
variability, inc l uding i ) thermonuclea r X-ra y bursts (e.g. 
Hoffman et "ail 1 19771 : IGallowav et~aH 120061 ; IZhang et all 
20 111 ), confirming the presence of a NS as the compact 



object, ii) millisecond oscillations during the X-ray bursts, 
indic a ting a spin freq u ency of ~581 Hz (|Zhang et al.l 
1 19971 ; IGiles et al. | l2002i ; IStrohmaver fe MarkwardtJ 120021 ) . 
and Hi) kHz quasi-periodic oscillations (kHz QPOs, e.g. 
Wiinands et al.1 1 19971; iBelloni et al.ll2007l ; lAltamirano et al.l 



20081 ; ISanna et alj |2012l ). From Eddington limited X-ray 



bursts, as suming a NS mass of 1 .4 M and a stellar radius 
of 10 km, IGallowav et al.l (|2006l ) estimated the distance to 
4U 1636-53 to be 6.0 ±0.5 kpc. Combining burst oscillations 
with phase-resolved optical spectroscopy of 4U 1636-53, 
ICasares et al.l (|2006l ) estimated the mass function and mass 
ratio to be f(M) = 0.76 ± O.47Af and q = 0.21-0.34, 
respectively, where q = M2/M1 with M2 the mass o f the 
donor and Mi the mass of the NS.[C asarcs et al.1 (120061 ) also 
showed that for a 0.48 Mq donor, the NS mass is 1.6-1.9 
Mq and the inclination i ~ 6 0-71°. The l atter result is 
in conflict with the model of iFrank et al.l (|l987h as the 
lack of X-ray dips in the light curve of 4U 1636-53 sets an 
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upper limit on the inclination of i < 60°. However, in this 
model the derived upper limit for the inclination depends 
on the location where the dips are formed, and since the 
exact geometry is not known the minimum inclination to 
observe dips could be slightly larger. On the other hand, 
the inclination could not be larger than ~75° since this 
would probably cause eclipses, which are not seen in this 
source. 

4U 1636-53 is a persistent X-ray source, although it 
shows variations in intensity up to a factor of 10, fol- 
lowing a narrow track in the colour-c olour diagram (CD 



and hardness-intensity diagram (HID; [Bclloni et al.l 



Altamirano et al ] |200SF ) with a ~40-d cycle ()Shih et al. 



2007; 



2005; 



Belloni et aT1l2007t ). The transition through the CD (or HID) 



is thought to be driven by changes in the mass accretion rate. 
The two distinct spectral states (hard, at the top right of 
the CD, and soft, at the bottom of the CD) are attributed 
to a different accretion flow configuration. In the soft state, 
where the X-ray intensity is high, the disc is hot and ionised, 
with an inner radius extending down to the ISCO (or down 
to t he stellar surface in case the ISCO is within the NS, see 
e.g. iDone et aljfeOQTh . As the intensity gradually decreases, 
the mass accretion rate eventually becomes too low to ionise 
the full accretion disc. The outer disc regions cool down, the 
viscosity decreases, and the overall mass accretion rate is re- 
duced. Eventually, mass accretion rate at the inner regions 
of the disc should decrease, leading to a recedin g inner edge 
and the X-ray intensity reaches its minimum |Done et al.l 
120071 ). The inner regions are likely replaced by a hot corona, 
and the spectrum hardens, which is the characteristic that 
defines the hard state. In the mean time, the secondary keeps 
transferring mass to the disc causing the density and tem- 
perature to increase, hence enhancing the viscous stress, and 
the mass accretion rate increases again. Eventually the disc 
gets ionised, the source rebrighten s and makes a transition 
to the soft state. IShih et al] l|2005h have shown that this re- 
brightening cycle for 4U 1636-53 of ~40 days is consistent 
with the viscous time-scale in the outer disc. 

High signal-to-noise and moderate-resolution spectra of 
4U 1636-53 r eveale d broad, asymmetric Fe emission lines. 
IPandel et all (|200Sl ) analysed three XMM-Newton/RXTE 
observations, once when the source was in the transitional 
state (between the hard and the soft state), and twice 
when it was in the soft state. In all three spectra they 
found that the Fe line profile is consistent with a relativis- 
tically broadened line coming from the inner disc, which 
appeared to be at the ISCO for a non-rotating N S. This 
finding is partly contested by ICackett et all (|2010l ) who fit 
the same spectra with a blurred reflection model that in- 
cludes effects like Com pton broadening and line emission 
from relevant elements. ICackett et al.l ( 2010T ) showed that 
fitting the spectrum with a self-consistent reflection model 
resulted in slightly larger values for the inner disc radius 
than when the sp ectrum was fit with the phe nomenological 
diskl ine model |Fabian et all 1 19891 ) u sed by | Pandel et al.l 
J2008I'). Nonetheless, both the results of IPandel et all (120081 ) 
and ICackett et all (|2010h are based on data which were not 
correc ted for pileup effects. As demonstrated by iNg et all 
(|2010h . re-analyses of the same XMM-Newton data consid- 
ering pileup and backgrou n d effe cts, suggested a different 
iron line profile. INg et all (|2010h found that the Fe lines 
in the three spectra of 4U 1636-53 used by IPandel et al.l 



|200Sft appeared to be symmetric and could be well fit with 
a Gaussian profile, although, due to the statistics, a rela- 
tivistic o ri gin o f the line could not be excluded. However, 
INg et all (|2010h did not include the simultaneous RXTE 
data in their fits, which can affect the Fe emission line profile 
that strongly depends on the underlying continuum. Trig- 
gered by the debate about the Fe line in 4U 1636-53, we 
obtained and analysed new XMM-Newton/RXTE observa- 
tions of 4U 1636-53, and re-analysed the previous XMM- 
Newton/RXTE data of 4U 1636-53. We investigated and 
corrected for instrumental effects, used a wide bandpass 
spectrum (0.8-120 keV), and fit the Fe emission line with 
the most commonly used phenomenological model and with 
a relativistically blurred ionised reflection model which self- 
consistently includes Compton broadening. Additionally, for 
the first time in Fe line studies, we investigated the contri- 
bution of different sources of disc illumination that could 
produce the reflection spectrum, and their contribution in 
the different spectral states. 

All together, 4U 1636-53 is an interesting source ex- 
hibiting several features that could be used to constrain key 
parameters like the NS mass and radius, as well as to test 
the standard NS accretion picture described by the colour- 
colour diagram. Comparing the properties of the Fe line as 
a function of spectral state - or even more specifically with 
the general properties of the spectrum like luminosity, black- 
body temperature of disc and boundary layer - may help us 
to constrain the origin of the line emission region. Further- 
more, kHz QPOs are seen to vary in freq uency depending on 
the spectral state of the source (see e.g. Belloni et al.l 1 20071 ; 
I Altamirano et alj|2008l ; ISanna etaHl2012h . Therefore, con- 
necting the position in the CD with the Fe line properties 
and the kHz-QPO frequencies, may help to break the degen- 
eracy in the models used to explain the spectral and timing 
features in accreting NS systems. In this paper, however, we 
focussed on the spectral analysis and the Fe line properties 
as this requires a careful investigation. In a companion pa- 
per (Sanna et al. 2013, in prep.) we combine the results of 
the spectral analysis presented in this paper with the timing 
analysis of the simultaneous RXTE observations. 



2 OBSERVATIONS AND DATA REDUCTION 

4U 1636-53 has been observed with XMM-Newton nine 
times in the last decade (between 2000 and 2009). The first 
two observations were taken with all CCD cameras oper- 
ated in imaging mode and were strongly affected by pileup. 
In the other seven observations, of all the CCD ca meras 
only the EPIC-pn (hereafter PN; IStriider et a.1.1 [200ll ') cam- 
era was on, operating in timing mode. In this mode one of 
the dimensions of the CCD is compressed into one single row 
to increase the read out speed. In our spectral analysis we 
only used the XMM-Newton data taken with PN in timing 
mode. We did not include the high-resolution data from the 
Reflection Grating Spectrometer in the energy range 0.33- 
2.5 keV. We note that the XMM-Newton observation taken 
on March 14, 2009 had a flaring high-energy background 
during the full ~40 ks exposure, and we therefore did not 
include this observation in our spectral analysis. 

Since 1996 RXTE observed 4U 1636-53 more than a 
1000 times; from March 2005 the source was regularly ob- 
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Table 1. XMM-Newton / RXTE observations of 4U 1636-53 used in our analysis 



Observation 


Instrument 


Identification Nr. 


Start Date 


Start Time 


Exposurc(ks) * 


Obs. 1 


XMM-Newton 


0303250201 


29-08-2005 


18:24:23 


25.7 




RXTE 


91027-01-01-000 




16:35:28 


26.2 (PCA) 
9.0 (HEXTB) 


Obs. 2 


XMM-Newton 


0500350301 


28-09-2007 


15:44:56 


14.3 




RXTE 


93091-01-01-000 




14:47:28 


26.9 (PCA) 
8.8 (HEXTB) 


Obs. 3 


XMM-Newton 


0500350401 


27-02-2008 


04:15:37 


34.7 




RXTE 


93091-01-02-000 




03:46:56 


25.3 (PCA) 
8.3 (HEXTE) 


Obs. 4 


XMM-Newton 


0606070201 


25-03-2009 


22:59:30 


23.8 




RXTE 


94310-01-02-03 




23:00:32 


1.9 (PCA) 




RXTE 


94310-01-02-04 


26-03-2009 


00:39:28 


1.6 (PCA) 




RXTE 


94310-01-02-05 




02:17:36 


1.4 (PCA) 




RXTE 


94310-01-02-02 




03:54:24 


1.3 (PCA) 
2.2 (HEXTE)" 


Obs. 5 


XMM-Newton 


0606070301 


05-09-2009 


01:57:03 


32.8 




RXTE 


94310-01-03-000 




01:17:36 


16.6 (PCA) 




RXTE 


94310-01-03-00 




08:20:32 


7.3 (PCA) 
7.6 (HEXTB)" 


Obs. 6 


XMM-Newton 


0606070401 


11-09-2009 


08:48:11 


21.1 




RXTE 


94310-01-04-00 




08:42:24 


18.4 (PCA) 



5.7 (HEXTE) a 

* Final exposure time after excluding X-ray bursts, detector drops, and background flares; see 
Section |2.1 1 and 1 2. 2 1 for more details. 

a Total exposure time of the combined HEXTE data; see Section |2. 2 1 for more information. 
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Figure 1. Long-term light curve of 4U 1636-53 based on pointed RXTE /PCA observations of the last ~5 years. The intensity is the 
2-16 keV Crab normalised PCA count rate (see Section l2.2.2l l. The numbered points mark the six XMM-Newton observations, with the 
numbers corresponding to the observations given in Table [T] 



served for ~2 ks every t wo days (for the firs t results of this 
monitoring campaign see iBelloni et aT]|2007l ). For our spec- 
tral analysis we only included the ten RXTE observations 
that were taken simultaneously with XMM-Newton. In Ta- 
ble [1] we give an overview of the observations, labelled in 
a chronological order. Using the pointed RXTE /PCA ob- 
servations of the last ~5 years, in Figure [T] we show the 
long-term light curve of 4U 1636-53 on which we marked 
the moments of the XMM-Newton observations used in this 
paper. 



2.1 XMM-Newton data reduction 

For the XMM-Newton data reduction we used the tasks of 
the Science Analysis System (SAS) version 12.1, applying 
the latest calibration files (xmmsas_20120523). To generate 
the concatenated and calibrated PN event files we used the 
meta-task epproc. We corrected for rate-dependent charge 
transfer inefficiency (CTI) effects using epfast, and we fil- 
tered out X-ray bursts, telemetry drops, and intervals of 
flaring high-energy background. For the latter, we excluded 
time-intervals in which the 10-12 keV count rate for single 
events (pattern=0) was larger than 1 cts s _1 . We gener- 



ated a time-averaged source spectrum per observation, us- 
ing single and double events only (pattern < 4), exclud- 
ing events at the edge of the CCD and at the edge of a 
bad pixel (flag=0), and using events within a 17-column 
wide box, centred on the source position, for which the cen- 
tral 3 columns were excised in order to correct for pileup 
effects (removing ~70% of the detected photons). We ex- 
cised the central 5 columns for Obs. 2 and Obs. 3 (remov- 
ing ~90% of the detected photons; see Appendix |X] for de- 
tails). We generated the response matrices using the task 
rmfgen. Following the recommendations given in the SAS 
User guide (Chapter 4.6.1 flfbr piled-up observations in tim- 
ing mode, we produced the ancillary response files using the 
task arfgen. For PN timing mode observations, the point 
spread function (PSF) of the telescope extends further than 
the CCD boundaries, and extracting the background spec- 
trum from the outer regions of the CC D leads to a spectrum 
contaminated with source photons (see lHiemstra et alj|201ll ; 
iNg et alll2010l . for a discussion on contaminated background 
spectra). To model the background, we extracted a spec- 
trum from the outer region of the CCD (RAWX in 5-20) 

1 http://xmm.esac.esa.int 
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of a "blank field" ob servation instead (see iNg et all l20ld: 
iHiemstra et al.|[201lf ). Based on similar sky coordinates and 
column density along the line of sight, we used the timing 
mode observation of the BH GX 339-4 (Obs 0085680601), 
in which the source was not significantly detected, as our 
"blank field" for all six XMM-Newton observations. Finally, 
we rebinned the source spectra such that we oversampled 
the PN energy resolution by a factor of 3 ensuring that we 
have a minimum of 25 counts per bin. 

2.2 RXTE data reduction 

2.2.1 Energy spectra 

For the RXTE data reduction we used the heasoft tools 
version 6.10, following the recipes in the RXTE cook boolQ. 
Applying the standard screening criteria and excluding time- 
intervals of detector drop-outs and X-ray bursts, we used the 
tool saextrct to extract the PCA spectra from Standard-2 
data, where we only included events from the third propor- 
tional counter unit (PCU2) being this the best-calibrated de- 
tector. We corrected for PCA deadtime, and applied a 0.6% 
systematic error to the PCA data. The PCA background 
was estimated with the tool pcabackest and the response 
files were generated using pcarsp. After excluding detec- 
tor drop-outs and X-ray bursts, we produced the HEXTE 
spectra for cluster-B events only, using the script hxtlcurv. 
With the tool hxtrsp we generated the HEXTE response 
files. No systematic errors were added to the HEXTE data. 
For the observations starting on March 25 and September 5, 
2009 (Obs. 4 and Obs. 5, respectively), RXTE did not cover 
the full XMM-Newton observation but several shorter expo- 
sures were taken instead. To speed up the fitting, we used 
the PCA spectrum with longest exposure, since all the spec- 
tra were consistent with each other within errors. However, 
in order to improve the poor HEXTE statistics we combined 
the individual HEXTE spectra with the tool sumpha. Net ex- 
posures of the final used PCA and HEXTE data are given 
in the last column in Table [T] 

2.2.2 Intensity and Colours of 4U 1636-53 

We calculated X-ray colours and intensity of 4U 1636-53 
using the Standard-2 PCA data. The intensity represents the 
2-16 keV count rate, and soft and hard colours are defined 
as the count rate ratio in the energy bands 3.5-6.0 keV / 
2.0-3.5 keV and 9.7-16.0 keV / 6.0-9.7 keV, respectively. To 
obtain the exact count rate in each of these energy bands 
we interpolated linearly in channel space. To correct for the 
gain changes and differences in the effective area between 
the PCUs, as well as to correct for the differences due to 
changes in the channel-to-energy conversion of the PCUs 
as a function of time, we normalised by the Crab Nebula 
val ues obtained close in time to our observations (for details 
see lKuulkers et al.lll99j ; lAltamirano et alj|2008h . We finally 
averaged the normalised colours and intensity per PCU for 
the full observation using all available PCUs. 

In Fig[2]we show the HID and CD of 4U 1636-53, with 
each point representing the average intensity, hard and soft 

2 http : //heasarc . gsf c .nasa . gov/docs/xte/recipes/cook_book 



colours per observation for a total of 1281 observations, up 
to May 2010. In these diagrams we also marked the position 
of the source at the time of the XMM-Newton observations. 
Based on the position of the source in both diagrams, we 
concluded that Obs. 1 and Obs. 6 were i n the hard state , 
also referred to as the transitional state (|Lin et al.l 120071 ). 
characterised by a relatively high hard colour and low X-ray 
intensity. The positions of Obs. 2-5 were consiste nt with 
the soft state (|van der Klisl 120061 ; iLin et all koOTT ). where 
both colours almost reach their mini mum values and the 
X-ray intensity is moderate. Following ILin et al.l (|2007l ). in 
the rest of this paper we refer to Obs. 1 and Obs. 6 as the 
transitional-state observations, and we refer to Obs. 2-5 as 
the soft-state observations. 

3 SPECTRAL ANALYSIS AND RESULTS 

For the spectral fitting of the XMM-Newton /RXTE data o f 
4U 1636-53, we used xspec version 12.7.1 (|Arnaudlll996T ). 
To study the evolution of the Fe emission line and its un- 
derlying continuum as a function of the spectral state, we 
fit each of the individual observations with a set of phe- 
nomenological models and a relativistically smeared reflec- 
tion model in addition to the direct emission components 
(see following Sections for the model description). In order 
to study the Fe line properly, we need to have a well-defined 
continuum on both sides of the Fe bandpass. For that rea- 
son we fit the XMM-Newton/RXTE spectra in the 0.8-120 
keV bandpass, with PN covering the 0.8-11 keV range (from 
channel 27 up to channel 267, unless otherwise mentioned), 
PCA taken in the 10-25 keV band (from channel 22 and 
53), and HEXTE covering the 20-120 keV range (between 
channels 11 and 45). Spectral uncertainties are given at 90% 
confidence (Ax 2 = 2.706 for a single parameter). 

In all our fits we included a phabs component to 
account for interstellar absorptio n, using the a bund ances 
and photo-el e ctric cross section of I Wilms et~ai1 (|2000l ) and 
IVerner et al.l (1 19961 ). respectively. Additionally, to account 
for flux calibration disparities between the different instru- 
ments, we added a multiplicative constant to the model. 

3.1 Spectral model 

3.1.1 Direct emission 

Spectral modelling of accreting NSs has been, and still is, a 
controversial matter. A variety of models have been pro- 
posed in the past years. All models, however, include at 
least a soft/thermal and a hard/Comptonized component 
that vary according with the source state. The soft /thermal 
emission comes from the accretion disk and from the NS 
surface, an d are typically model led with a multicolour disk 
blackbody (|Mitsuda et al.l 1 19841 ) and a blackbody, respec- 
tively. The hard/Comptonized emission is probably created 
by the up-scattering of soft photons coming either from the 
accretion disk, the NS surface, or both simultaneously, in 
the hot electron gas surrounding the central region of the 
system (usually referred as corona). The hard emission is 
modelled with a power law with or without cut off at high 
energies, although thermal comptonisation models are also 
.html widely used. 
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Figure 2. Hardness-Intensity diagram (left) and colour-colour diagram (right) of 4U 1636-53. Each grey point represents the averaged 
Crab- normalised colours and intensity per observation, for a total of 1281 observations (see text in Section 12.2.21 for details). The 
numbered black points mark the position of the six XMM-Newton observations, with the numbers corresponding to the observations 
given in Table [TJ 



In this work we used a multicolour disc blackbody 
(diskbb) to fit the thermal emission from the accretion 
disc, which was required by the fits, both in the soft and 
the transitional state. To model the thermal emission from 
the NS surface we used a single-temperature blackbody (see 
iLin et al]|2007f ). For the hard emission we used the thermally 
comptonized continu um model nthcomp (|Zvcki et al.ll 1999 ; 
IZdziarski et al.lll996l '). which compared to an exponentially 
cutoff power law describes more accurately the high energy 
shape and the low energy rollover, with similar number of 
parameters. This model allows to select the spectral shape 
of the source of seed photons between a (quasi)blackbody 
and a disk blackbody. During the fitting process we used 
either of the thermal components as source of seed pho- 
tons for nthcomp. Both options gave statistically accept- 
able fits, however when we linked the emission from the NS 
surface/boundary layer to nthcomp the blackbody emis- 
sion itself turned out to be negligible. We therefore opted 
for the disk blackbody as the source of seed photons for the 
nthcomp component. 

Additionally, the emission coming from the inner parts 
of the accretion disc, close to the NS, may be modified by 
relativistic effects. For a rotating NS, thes e effects are similar 
as for a spinning BH (|Miller et al.l Il99^ ) and are approxi- 
mately described by a Kerr metric characterised by the di- 
mensionless angular momentum, a, = cJ/(GM 2 ). For the 
case of a NS, the metric also depends on the internal struc- 
ture of the NS, set by the equation of state (EoS), which 
is unknown. However, the spin par ameter for a NS c an be 
approximated as a* = 0.47/Pf msl dBraie et alj|2000h. For 
4U 1 636-53, with v = 581 Hz |Strohmaver fc MarkwardrJ 
120021 ). the spin paramet e r is ^ 0.27. This value is consistent 
with what iMiller et all \ 19981 ) found where, for a 500 Hz 
rotating NS, a given EoS, and a NS mass of 1.4 M©, the 
spin parameter is ~0.23. Scaling this value up to a rotation 
of 581 Hz gives a, ~ 0.27. To test whether the direct disc 
emission in the spectrum of 4U 1636-53 is significantly af- 
fected by relativistic effects we tried, instead of the diskbb, 
a multi-temperature blackbody model for a thin, general rel- 
ativistic accretion disc in a Kerr metric, which also includes 



self-irradiation (kerrbb; iLi et aT]|2005l ). However, using the 
kerrbb, assuming a distance of 6.5 kpc, a NS mass of 1.4 
M©, and a* = 0.27, did not significantly improve the fit. 
Moreover, the uncertainty in the distance and the NS mass 
of 4U 1636-53, makes the outcome of this model less reli- 
able. Therefore, for the rest of our analysis we only used the 
diskbb component to fit the direct disc emission. 

To summarise: we fitted both the soft and the transi- 
tional state of the source with a multicolour disk blackbody 
(diskbb), plus a single-temperature blackbody (bbody) and 
the thermally comptonized component (nthcomp) with the 
soft seed photons coming from the accretion disk. 



3.1.2 Phenomenological and reflection models for the iron 
emission line 

Phenomenological models: 

After fitting the continuum we found several residuals in 
the whole PN range, with the most prominent one in the 
energy range 4-9 keV around the Fe line emission region. 
In order to fit these residuals we added a Gaussian emis- 
sion line with the energy constrained between 6.4 keV and 
6.97 keV (Fe K band), while the other parameters of the 
line (normalisation and width) were free to vary. In all six 
observations the data were well fitted by this component, 
however all the lines were very broad, showing widths rang- 
ing from ~1.0 to ~1.4 keV (see Table [6] for more details on 
the fit parameters) . Triggered by this an d by previous results 
on the Fe emission line in 4U 1636 -53 (|Pandel et al.ll2008l ; 
ICackett et"aHl20ld ; iNg et al.l [20T0I ) . we used a set of phe- 
nomenological models describing relativistically-broadened 
lines. We selected three of the most commonly used models 
by the community: diskline, laor, and kyrline. 

disklin e is a relativistic mo del for a Schwarzschild met- 
ric, a* = 0, (|Fabian et all 1 19891 ). thought to be suitable for 
NS with dimensionless angular parameter lower than 0.3 
for which the met ric should only m arginally deviate from 
Schwarzschild (see lMiller et al.lll998l . for details on the sub- 
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ject). Notice that the diskline model does not include light 
bending corrections. 

LAOR is a relativistic line model which includes the 
effects of a maxim ally rotating Kerr black hole, a* = 
0.998, (|Laorl Il99ll ). including light bending corrections. 
Although this model was meant for black holes, it has 
been largely used to fit disc lin es in NS systems (see 
iBhattacharvva fc Strohmaverll20"07l . for comparison between 
diskline and LAOR in the NS LMXB Serpens X-l). 

The kyrline co mponent includes all relativistic effects 
for the Kerr metric jDovciak et al.l [20041 '). Differently from 
the previous two models, the metric changes with the spin 
parameter, and can be adjusted for any value of a* between 
and 0.998. The line profile is calculated for emission from 
outside the ISCO only, and includes the effects of limb 
darkening. 



Reflection model: 

When X-rays irradiate an optically thick material such as 
the accretion disc, they produce a reflection spectrum in- 
cluding fluorescence lines, recombination and oth er emission 
(see e.g. lFabian et al.ll 19891 ; iRoss fc Fabianll 19931 ). When the 
illuminating flux is high enough (possibly combined with the 
thermal blackbody radiation intrinsic to the accretion flow) 
the surface of the accretion disc is expected to be ionised. 
The ionisation state of the reflecting material determines th e 
shape of the reflection spectrum (e.g. IRoss fc Fabianll 1993T ). 
and thus it is important to solve for the thermal balance of 
the disc. In most X-ray sources the incident emission for the 
reflection spectrum is generally a hard power-law spectrum. 
However, in NS systems, the emission coming from the NS 
surface/boundary layer may contribute as well. As the shape 
of the reflection continuum also depends on the incident 
emission, it is important to investigate the different pos- 
sible sources of irradiation. Several ionised reflection models 
are publicly available within xspec, which self-consistently 
compute the ionisation and thermal balance according to 
the radiation field. In our analysis we used simultaneously 
two different reflection models: rfxconv and bbrefl. 

rfxconv is an u pdated ver sion of the code in 
iDone fc Gierlinskil |200rj ). using the IRoss fc Fabianl (120051 ) 
reflection atables. This model can be used with any input 
continuum and has therefore the advantage of not having a 
fixed exponential cut-off in the illuminat i ng po wer law at 300 
keV as in the model of IRoss fc Fabianl l|2005l ). We decided 
not to reflect the blackbody emission with rfxconv because 
this component calculates the ionisation balance assuming 
that the ionising spectrum is a power law, and this is not ap- 
propriate for the case of the reflection of blackbody emission 
off the disc. Instead, we used bbrefl, which provides the re- 
flection spectrum from a constant density disc illuminated 
by a blackbody (jBallantynd 12004] ) . Both components also 
include lines and edges of the most important elements and 
ions which are Compton and thermal broadened depending 
on the ionisation state of the reflecting material. 

Since the reflection spectrum may be further smeared 
by relativistic effects in the inner regions of the accre- 
tion disc, we convolved the reflected emission with a Kerr- 
metric kernel which includes r elativistic effects (kerrconv; 
iBrenneman fc Reyn olds 2006), where we fixed the spin pa- 



rameter to 0.27. In this work we explored how the corona 
and the boundary layer contributed individually to the re- 
flection spectrum. Details on the specific parameters and 
their settings are given in Section 13.21 and the results are 
discussed in Section [3. 41 



3.1.3 Residual absorption and emission features 

After fitting the spectra with the continuum model plus 
the relativistic Fe line profile, the XMM-Newton data still 
showed some narrow residuals, especially at low energies. 
Besides the commonly seen residuals at the instrumental Si- 
K and Au-M edges at ~1.8 keV and ~2.2 keV, respectively, 
there were a few other clear features apparent in the residu- 
als. All six observations showed absorption features at ~0.9 
keV, at ~1.5 keV and at ~ 9.2 keV. Taking into account 
the moderate PN energy resolution (~70 eV at 1 keV), the 
~0.9 keV feature could be due to absorption by Ne IX (0.905 
keV) or Fe XVIII (0.873 keV), and may possibly be of astro- 
physical origin. The energy of the ~1.5 keV feature is either 
consistent with the Fe XXI-XXIV blend or with the Al-K 
edge. The effective area curves of the PN show a strong Al-K 
edge at ~1.56 keV, which suggests that this feature is likely 
an artefact of the PN calibration. The absorption feature at 
~9.2 keV is consistent with the Fe XXVI edge. The above 
residuals contributed to an enhanced \ 2 - To improve the fit 
we modelled these residuals with a Gaussian absorption or 
emission component. Since the main focus of this paper is 
the spectrum at energies near and above the energy of the Fe 
line, we do not discuss these lines further. In Obs. 5 we fur- 
ther ignored 11 channels from the PN data (spread between 
channels 50 and 148, corresponding to the energy range ~ 
1.3 to ~ 4.7 keV), that appeared in the spectrum like very 
narrow absorption or emission features. 



3.1.4 Fitting procedure 

We started fitting the energy range 0.8-100 keV using PN 
(0.8-10 keV), PCA (3-20 keV) and HEXTE (20-100 keV) 
data. From the fits with the phenomenological model kyr- 
line for the line (see Sec. 3.1.1 for details about the model) 
we noticed a mismatch between PN and PCA data between 
4 and 8 keV. To investigate the mismatch we allowed the 
phenomenological model of the line to vary between PN and 
PCA, and we found that, except for the normalisation (with 
the iron line being a factor ~ 4 stronger in the PN than in the 
PCA spectra), the parameters of the Fe line were consistent 
within errors. By letting the line normalisation free between 
the two instruments the fit improved significantly. The val- 
ues of the normalisation in PN and PCA were, in both cases, 
significantly different from zero. We tested the fits excluding 
both PCA and HEXTE data, and we fitted the model pre- 
viously described to the PN data only. We found that, for 
the observations in the transitional state the line parameters 
do not change much between PN, and PN+PCA+HEXTE, 
except for the normalisation of the line and the inner radius 
of the disc that are marginally larger in the case we used 
the three instruments. All continuum components were sta- 
tistically required to fit the PN-only spectrum. In the soft 
spectra all the line parameters were consistent within the 
errors with being the same. In this case the direct emission 
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from the corona is not statistically required by the fit. This 
is however understandable given that the PN only extends 
up to 10 keV and the hard emission is less important in the 
soft than in the transitional state. 

Based on these comparisons, and to be able to properly 
constrain the underlying continuum and model the Fe line 
emission, we proceeded further by including all three instru- 
ments in our fits, excluding the PCA data below 10 keV, 
where we relied on the PN data. We further modelled the 
residual absorption between 9 and 10 keV using an edge with 
energy free between 9.2 and 9.3 keV, which could be due to 
Fe XXVI. The optical depth of the edge varied between 0.05 
and 0.1. Not adding the edge to the fits does not affect the 
parameters of the line. The high level of ionisation required 
to create this feature contrasts with the Fe line energy val- 
ues found from the fits, suggesting that this absorption edge 
might be due to a calibration mismatch between PN and 
PCA. 



3.2 Parameter settings 

3.2.1 Continuum, 

As we mentioned, we used the diskbb, bbody and nthcomp 
components to model the direct emission coming from the 
accretion disc, boundary layer, and corona, respectively. The 
parameters for these components are the disc temperature 
at the inner edge, fcTm, the temperature of NS surface, kTsB , 
and the Comptonisation photon index, the electron temper- 
ature (high energy rollover), and the seed photon tempera- 
ture (low energy rollover), F, kT e , and fcTdbb, respectively. 
A^dbb, Nbb, and -/Vnth are the normalisation parameters for 
the diskbb, bbdoy and nthcomp components, respectively. 
All the direct emission parameters were free to vary (except 
fcTdbb which was set equal to kT[ n of diskbb), but were cou- 
pled between the different instruments. 

3.2.2 Phenomenological line models 

The phenomenological line models are parameterised by the 
inclination angle of the accretion disc, i, the rest energy of 
the line, Eu ne , the inner and outer edge of the disc, 7?i n and 
Rout, the emissivity index, /3 (kyrline has the option to 
allow the emissivity index to be different in the inner and 
outer disc regions; here we chose to have a disc described by 
a single emissivity index) , and the normalisation of the line, 
Mine, in photons cm -2 s _1 . We constrained E\- me to range 
between 6.4 and 6.97 keV, and fixed the outer disc radius 
to be 1000 GM/c 2 . Using kyrline we fitted the data fixing 
a* to three different values: 0, 0.27 and 0.998, and limited 
the inner radius to the ISCO R in > 6, 5.12, 1.23 GM/c 2 , 
respectively. 



3.2.3 Reflection model 

The ionised reflection emission is characterised by three pa- 
rameters: the scaling reflection factor, Sl ro fl, from rfxconv 
(f2 re fl < represents only the reflected component), bbrefl 
normalisation (flux per emitting area), and the ionisation 
parameter, £ = 47tF/tih , with F the total illuminating flux, 
and riH the hydrogen number density. Under the assump- 
tion that the reflection region in the accretion disc has a 




Figure 3. The two transitional-state spectra (top, Obs. 1, and 
bottom, Obs. 6) of 4U 1636-53. Each plot shows the simultane- 
ously fitted XMM-Newton/RXTE spectra and unfolded model 
in the main panel, and the residuals in units of sigmas in the 
sub-panel. To account for the direct emission from the accre- 
tion disc, NS surface/boundary layer, and corona, this model 
contains a diskbb (black/line), BBODY (blue/dashed), and NTH- 
COMP (black/dashed-dottcd) component, respectively. The rela- 
tivistically smeared reflected emission is due to two incident com- 
ponents: the NS surface/boundary layer (cyan/dotted) and the 
corona (light grey/dashed-triple dotted). 



constant ionisation parameter we coupled the ionisation pa- 
rameter between rfxconv and bbrefl. 
The relativistic effects, determined by the metric and the 
disc properties, were included using the kerrconv compo- 
nent which is parameterised by the spin parameter, a,, the 
disc inclination, 0, the disc inner and outer radius, Ri n and 
-Rout (in units of the marginal stable radius, i? ms ), and the 
disc emissivity index, f3. The latter could be different for the 
inner and outer disc, although we chose to have a disc de- 
scribed by a single emissivity index. Ri n , 6, and j3 were free 
to vary, but we fixed R ou t to its default value (400 i? ms ) and 
a, to 0.27 (see Section |3.1.1| I. G iven t his spin parameter, us- 
ing equation (3) in lMiller et all (|l998T l, the ISCO is at ~5.T2 
gravitational radii, R g (R g = GM/c 2 , this is the first-order 
approximation valid for spin frequencies below ~400 Hz). 
All free parameters in the relativistically smeared reflection 
model were coupled between the different instruments. 
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The inclination in rfxconv is given as a cosine function 
and we fixed cos to a value consistent with the inclination 
angle in kerrconv obtained from the initial fits of the rel- 
ativistically smeared reflection model. For all observations 
this initial value for the inclination was in the range 63-69° , 
which is quite high but still consistent with the results of 
ICasares et al.l (|2006T ). However, for most of the observations 
(except Obs. 2 and 4) we found that the best-fitting values 
for the inclination inferred from kerrconv were larger than 
80°, which is too high as no eclipses have been seen in this 
source. In Section 14. II we examine this issue in more detail. 
Besides the parameters we directly got from the blurred re- 
flection models, it is also common to provide the strength 
of the Fe line, usually expressed in terms of its equivalent 
width (EW). However, as the Fe line is part of the reflec- 
tion spectrum, it is not possible to separate the emission line 
from the underlying reflection continuum. Therefore, we es- 
timated the equivalent width of the reflection continuum in 
the region of the Fe line. As the line is broadened, we used 
the 4-9 keV range as our line region for all six observations. 
The equivalent width is defined as flux/contin; with flux be- 
ing the (unabsorbed) 4-9 keV flux of the reflection spectrum, 
and contin the average (unabsorbed) 4-9 keV flux density 
of the direct continuum. In Table [2] we give the equivalent 
width obtained in this way. However, we note that these val- 
ues probably overestimate the true equivalent width of the 
Fe line since flux not only includes the line photons but also 
the photons of the underlying reflection continuum. 

3.3 Sources of irradiation 

To investigate how the boundary layer and the corona con- 
tributed to the reflection spectrum in 4U 1636-53, we de- 
fined our model as: phabs*(diskbb+ bbody + nthcomp 
+ kerrconv* (bbrefl + rfxconv*nthcomp) ) (hereafter, 
model 1), where both the bbody and nthcomp components 
were reflected and relativistically smeared. We found accept- 
able fits for the observations, with the reduced \ 2 ranging 
from 1.06 to 1.12 for 278 d.o.f. However, for some observa- 
tions, we found that one of the two reflection components 
was not contributing significantly to the total spectrum. We 
therefore investigated the scenario where only one of the two 
incident components, bbody or nthcomp (hereafter model 
2 and 3, respectively), contributed to the reflected emission. 

In Table [3]we summarised the values of x 2 an d number 
of degrees of freedom for the different fitting models applied. 
For three out of the 6 observations (Obs. 1, 3 and 5) the best 
fit was obtained by including both reflection from the NS 
surface/boundary layer and the corona (model 1). For Obs. 
2 and 4, both model 1 and 2 showed the same \ 2 with only 
1 d.o.f. of difference. f2 re fl, in Obs. 2 and 4, was consistent 
with being zero, hence rfxconv*nthcomp in model 1 was 
not required to model the data. For Obs. 4 we then decided 
to use model 2. For Obs. 2, on the other hand, we finally 
opted for model 3. Fitting results from model 2 showed that 
the direct emission from the NS surface/boundary layer was 
consistent with zero within errors. This would imply that 
the NS surface/boundary layer, although not directly visi- 
ble, still contributes to the reflection spectrum. While this 
is in principle possible, it is more likely that the model used 
to fit this spectrum has too many components, resulting 
in an over-parameterisation of the problem. Taking all this 



into account, for Obs. 2 we discarded models 1 and 2, al- 
though statistically preferred, and we selected model 3 in- 
stead. The case of Obs. 6, was very similar to that of Obs. 2 
(see [3]l . Also in this case, fitting model 1 we found that 
only the bbrefl component was significantly contributing 
to the reflected emission, while switching to model 2 the 
direct emission from the NS surface/ boundary layer was 
not significantly required. As for Obs. 2, we decided to use 
model 3. In Figure[3]we show the spectra and the individual 
components of the best-fit models for the transitional obser- 
vations fitted with the reflection model. The top and bottom 
panels of Figure[3]show Obs. 1 and 6, respectively. As previ- 
ously mentioned, in Obs. 1 both the reflection from the NS 
surface/boundary layer (cyan-dotted line) and corona (light 
grey dashed-triple dotted line) contribute significantly, with 
the latter accounting for ~30% of the reflection emission in 
the 4.0-9.0 keV range. On the other hand, in Obs. 6 the 
NS surface contribution to the reflection spectrum is not re- 
quired. 

In Figure [4] we show spectra and individual components of 
the best-fit models for the soft-state observations (top-left, 
Obs. 2, top-right, Obs. 3, bottom-left, Obs. 4, and bottom- 
right, Obs. 5). Obs. 2 is the only soft-state observation 
where the reflected spectrum is entirely described by reflec- 
tion from the corona (light grey dashed-triple dotted line). 
Obs. 3, as well as Obs. 5, show contribution from both re- 
flection components, with the the corona (light grey dashed- 
triple dotted line) contributing ~15% and ~13%, respec- 
tively in the Fe line region. Finally, Obs. 4 is the only soft- 
state observation where the corona does not contribute to 
the reflected spectrum. 



3.4 Evolution of continuum parameters and line 

3.4-1 Continuum parameters 

The continuum parameters from the reflection and phe- 
nomenological models are shown in Table [2] and Tables [4][5] 
respectively. In order to study the evolution of the individ- 
ual spectral components in relation to the position of the 
source in the CD, we ordered the observations according to 
their colour-colour coordinate S z , which is t hought to be a 
funct ion of the mass accretion rate (see e.g.. iKuulkers et al.l 
Il994h . The order of the six XMM-Newton/RXTE observa- 
tions according with their 5* z value is: 1-6-2-3-5-4. Apply- 
ing this order, in Figure [S] we show how several of the spec- 
tral parameters changed with the position in the CD, con- 
sidering both the reflection model and a simple relativistic 
emission model (kyrline with angular parameter fixed to 
0.27). From the first panel of Figure [5] we found that the 
disc temperature increased as the source became brighter, 
consistent with t he standard accretion disc model (see e.g., 
iDone et al.ll2007l . and reference therein). The disc flux (sec- 
ond panel) did not show a similar trend as the disc tem- 
perature, however it is apparent that the flux is on average 
lower in the transitional-state (Obs. 1 and 6) than in the 
soft state; although that is more evident when using the 
phenomenological model (red stars), it can still be noticed 
for the reflection model. From panels 3 and 4, we found that 
the flux and temperature of the NS surface/boundary layer 
did not show any clear trend with 5 Z in the fits with the 
reflection model. Notice, however, that the blackbody flux 
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Table 2. Best-fitting results for the reflection model. 
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NOTES.— A * denotes that the error has reached the maximum or minimum allowed value of the parameter. Uncertainties are given 
at a 90% confidence level. f! re fl < represents reflection only, where the actual reflection normalisation is |O rc fl|. N^bb ls defined as 
{Rin/ D\a) 2 cosd, with Ri n in km, D10 the distance in 10 kpc, and the inclination angle of the disc. JVgB is L3g/D 2 , where L%g is the 
luminosity in units of 10 39 erg s — 1 . %th is in units of photons keV - 1 cm -2 s — 1 at 1 keV. All the flux values reported represent the 
unabsorbed flux in the energy range 0.5-130 keV, all in erg cm -2 s —1 . EW represents the equivalent width of the reflection component in 
the 4-9 keV energy range. 



in the phenomenological model showed a clear increase be- 
tween the two source states. The bbody behaviour across 
the CD does not seem to follow the standard scenario, for 
which an increase in mass accretion rate should heat up the 
surface of the NS. Panels 5 to 7 show the corona emission 
properties across the CD. Both, photon index and flux of the 
component did not vary systematically as the source moved 
from transitional to soft state. Notice, however, F from the 
reflection model in average larger in the soft state. On the 
other hand, the electron temperature (kT e ) of the scatter- 
ing cloud clearly decreased as mass accretion rate increased 
(with the exception of Obs. 4 in the reflection model). The 
latter result is consistent with the scenario in which, as the 
source moves from transitional to soft state, the accretion 
disc becomes hotter emitting a larger number of soft photons 
that will cool down the electrons in the corona. 



3.4-2 Emission line: phenomenological models 

The line parameters from the phenomenological models are 
given in Table [75] In Figure [7] we show three representative 
spectra and individual components of a transitional-state 
observation {top, Obs. 1), and two soft-state observations 
{middle, Obs. 3, and bottom, Obs. 4) fitted with relativistic 
model kyrline with a* = 0.27. 

The first thing to notice from the table is the fact that 
the Fe line is very well modelled with a simple, symmetric, 



gaussian profile. Replacing the gaussian profile with a rel- 
ativistic line model always lead to an increase of \ 2 (with 
the exception of Obs. 4), with a maximum Ax 2 between the 
gaussian and the other models of ~32 (for 2 d.o.f. difference). 
The gaussian component, however, showed always a very 
broad profile, with a line width between 1 and 1.4 keV that 
seemed to increase as the source went from the transitional- 
state to the soft state, with the exception of Obs. 4 in which 
the width dropped to its lowest value of 1 keV. In Obs. 1, 
2, 5 and 6, all the model consistently showed an energy line 
at 6.4 keV that implies neutral or lowly ionised iron. In 
Obs. 3, kyrline gave energy values of ~6.7 keV while the 
other models place the line at ~6.4 keV. Finally, in Obs. 4 
diskline and kyrline gave an energy value of ~6.4 keV, 
while GAUSSIAN and LAOR gave values larger than 6.7 keV. 
The inclination was generally too high given that no dips or 
eclipses have been observed in 4U 1636-53. All the models 
showed inclinations larger than 70 degree, with most of the 
values showing an upper confidence limit that pegged at 90° . 
We further found that the emissivity index of the disc was 
roughly the same for all observations for all the relativistic 
models. 

In Figure [fJJ from top to bottom, we show the evolution of 
the EW, inner radius and line flux across the CD for the 
different components used to model the Fe K-q line. GAUS- 
SIAN and LAOR showed the largest EW, although with big 
errors. Within errors, we found the values of the EW to 
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Figure 4. The four soft-state spectra (top-left, Obs. 2, bottom-left, Obs. 4, top-right, Obs. 3, and bottom-right, Obs. 5) of 4U 1636-53. 
Each plot shows the simultaneously fitted XMM-Newton/RXTE spectra and unfolded model in the main panel, and the residuals in 
sigma units in the sub panel 

. Lines and colours are the same as in Figure [3] 



be consistent between the different models, and not vary- 
ing with the source state. In the middle panel of Figure [6] 
we show the inner radius in units of the gravitational ra- 
dius, R g . In the panel we plot 3 radii representing the ISCO 
for the Schwarzschild metric (6 Rg), the Kerr metric with 
a, — 0.27 (~5.1 R g ) and the Kerr metric with a, = 1 (~1.2 
R g )- Interesting to notice is that only in Obs. 1 all the phe- 
nomenological models gave a consistent estimate of the in- 
ner radius of ~ 11 R g , while for the other observations the 
different models showed significantly different inner radii. 
Among the models, LAOR showed generally the lowest val- 
ues of the radius. Although LAOR and kyrline with a* = 1 
are based on the same Kerr-metric kernel, in 4 out of 6 ob- 
servations both gave significantly different values of Ri„. A 
similar consideration is valid also for diskline and kyrline 
with a, = 0. Overall, the inner radius did not vary much 
from Obs. 6 to Obs. 3, and than significantly dropped in the 
last two observations. 

Finally, the bottom panel of Figure [5] shows the unabsorbed 
line flux in the 0.5-130 keV energy range. The GAUSSIAN 
and LAOR components showed the largest flux values in all 
observations. Furthermore, the flux values of both these line 
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BBREFL 


RFXCONV 
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313/278 


326/279 


334/279 
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286/278 


286/279 


288/279 
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285/278 
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301/278 


301/279 


327/279 
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287/268 


302/269 


308/269 


6 


290/278 


290/279 
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Table 3. \ and number of degrees of freedom for 3 different 
configurations for the incident emission. BBREFL models a reflec- 
tion spectrum from a disc illuminated by a blackbody. RFXCONV 
is used to create reflection assuming NTHCOMP as incident spec- 
trum. The third option includes both contributions from BBODY 
and NTHCOMP direct emissions. In grey we highlighted the best fit 
models chosen for the analysis. See Section l3.2.3l for more details. 



models increased from Obs. 6 to Obs. 3, and then decreased 
down to flux values comparable to the transitional-state 
ones, diskline and kyrline showed a more clear trend with 
St,: the line fluxes increased from the transitional-state to 
the soft state with the only exception of Obs. 4 where the 
flux dropped. 
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Figure 5. Evolution of several spectral parameters as the source 
moves through the CD. The observations are ordered according to 
their S z coordinate. From top to bottom, the panels show respec- 
tively the evolution of the DISKBB temperature, the disc flux, the 
BBODY temperature, the blackbody flux, NTHCOMP photon index, 
NTHCOMP flux and NTHCOMP electron temperature. The fluxes are 
unabsorbed in the 0.5-130 keV range in erg cm" 2 s~ 1 . 



3-4-3 Emission line: phenomenological models vs. 
reflection 

As shown in section 3.4.2, the iron emission line is well fitted 
by a symmetric gaussian profile, characterised however by a 
large breadth. In order to explore the broadening mechanism 
that shapes the reflection spectrum, we started by fitting 
the reflection spectrum without the kerrconv kernel that 
models the relativistic effects. The model therefore was (see 
Sec. 13. 3ft : phabs*(diskbb + BBODY + NTHCOMP + bbrefl 
+ rfxconv*nthcomp). The only broadening mechanism in 
this case is Compton broadening. We found that in all the 
six observations the model did not fit well the data, with 
xt varying between 1.3 and 1.7 (with d.o.f. ranging from 
271 to 282). From the residuals was clear a strong wiggled- 
like feature in the energy range 5 to 8 keV, consistent with 
the unfitted iron emission line. By adding the kerrconv 
component all the fits improved, as reported in Table [2] We 



found that the inner radii inferred from the reflection model 
are consistent with the values from kyrline with a* = 0.27. 
Furthermore, the reflection model gave very large inclination 
values, generally larger than 80°, consistent with kyrline 
(with the only exception of Obs. 2 in which the reflection 
model gave an inclination of ~50°). Finally, we noticed that 
the EW in the 4-9 keV energy range from the reflection 
model was consistent within errors with the EW from the 
relativistic model kyrline. 



4 DISCUSSION 

We fitted the X-ray spectra of six simultaneous RXTE and 
XMM-Newton observations of the NS low-mass X-ray binary 
4U 1636-53. For the first time we used a reflection model in 
which simultaneously both the NS surface/boundary layer 
and the corona were taken to be the source of irradiation to 
investigate the relative reflection of each of them. We further 
modelled the Fe emission line assuming different line profiles, 
both symmetric (GAUSSIAN) and relativistically-broadened 
(diskline, laor, and kyrline). 

We found that, both in the transitional state and the 
soft state, the corona is the most prominent component 
in the broadband spectrum. The emission from the NS 
surface/boundary layer is on average contributing an or- 
der of magnitude less than the corona to the continuum 
emission. However, the relative contribution of the NS sur- 
face/boundary layer increases going from the transitional to 
the soft state. In four out of six observations (Obs. 1 in the 
transitional state, and Obs. 3, 4 and 5 in the soft state) 
the NS surface/boundary layer is the main source of disc 
irradiation, accounting for ~70% and ~80% of the reflected 
spectrum in the Fe line region (based on the 4-9 keV unab- 
sorbed flux) in the transitional and soft state, respectively. 
The NS surface/boundary layer reflection spectrum is even 
more predominant in Obs. 4, where the corona does not 
contribute at all to the reflection. This may be due to the 
NS surface/boundary layer being the main contributor to 
the reflection spectrum simply because is closer to the inner 
edge of the accretion disc where most of the disc emission 
generates. On the other hand is also likely that part of the 
photons emitted by the NS surface/boundary layer interact 
with the corona by means of Compton scattering, as a con- 
sequence of that what we model as direct emission from the 
NS is only an underestimation. 

For the other two observations (Obs. 2 in the soft state, 
and Obs. 6 in the transitional state) the reflec tion spectrum 
origin ates entirely from the corona emission. ICackett et alj 
(|2010T ) analysed Obs. 1-3 using a blurred reflection model 
where only the NS surface/boundary layer (fitted with a 
blackbody component) was assumed to illuminate the disc. 
We found that the line properties inferred from our fits 
with the refle ct ion m odel are consistent within errors with 
ICackett et al.l (|2010n findings. The only discrepancy be- 
tween the results is related to the logarithm of the ioni- 
sation parame t er for Obs. 1, which we find to be ~1 while 
ICackett et al. (2010) find an ionisation parameter of ~2.5. 
It is not surp rising, at least for Ob s. 1 and 3, that our re- 
sults matched ICackett et al.l (|2010T l findings. As mentioned 
above, in Obs. 1 and 3, the NS surface/boundary layer is 
the predominant source of disc irradiation in the Fe line re- 
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Figure 6. Evolution of the Fe emission parameters of different line models as the source moves through the CD. The observations are 
ordered according to their 5 Z coordinate. From top to bottom, the panels show respectively the evolution of the equivalent width, the 
inner radius in units of the gravitational radius, R g = GM/c 2 , and the line unabsorbed flux in units of 10 — 11 erg cm -2 s _1 . Dashed, 
dashed-dot, and dashed-dot-dot lines (middle panel) represent the 6 R g , 5.1 R g , and 1.2 R s radii, respectively. 



gion, accounting for 70% and 85% of the reflected spectrum, 
respectively. On the other hand, it is interesting to mention 
that although in Obs. 2 we modelled the reflection spec- 
trum assuming the corona as source of irradi ation, the line 
param eters are still consistent with those of ICackett et al.l 
(|20ld) . 

We modelled the Fe emission line assuming different 
line profiles, both symmetric (GAUSSIAN) and relativistically- 
broadened (diskline, laor, and kyrline). We found that 
all the models fit reasonably well the data (see Table[B]for all 
the fitting parameters). The GAUSSIAN profile seemed to fit 
the line better, which suggests a sym metric prof i le for the Fe 
emission line. This result agrees with lNg et all (j201fjh . that 
found the lines from Obs. 1-3 to be symmetric, and that the 
GAUSSIAN and the relat i vistic profile LAOR could fit equally 
well the line. iNg et alJ (|2010l ) suggested that, contrary to 
previous claims, the width of the line can be explained by 
mechanisms oth er than relativistic effect s, such as Comp- 
ton broadening. iKallman fe White! (jl989h showed that iron 
lines generated in the ADC are characterised by a symmetric 
profile due to blending, Compton scattering and rotation. 
IKallman fe White! \ 19891 1 further estimated that Fe lines 
generated in the ADC could show EW values up to 100 eV 
assuming standard parameters. From the fits with the GAUS- 
SIAN model we found EW values betwe en ~300 and ~700 
eV, which are significantly larger than IKallman fe White] 



(jl98fj! ) estimates. This indicates that Compton scattering 
alone cannot explain the EW from the GAUSSIAN model. 
This is also confirmed by the reflection model: We find that 
excluding the convolution component kerrconv, the re- 
flection model fails to fit the Fe line in the spectrum. For 
complete n ess w e should mention that the EWs reported by 
INg et all (|2010f ) are significantly lower (between ~30 and 
~200 eV) than the ones we report here, and that INg et alJ 
(|2010l ) raised doubts on the Fe line with the largest EW 
(Obs. 1) being realistic. The fact that we use a wide band- 
pass energy spectrum (0.8-120 keV) to constrain the direct 
emissio n explains the d ifferences between our results and 
those of lNg et all <j2010T ). That said, our results suggest that 
even though the Fe emission line profile is compatible with 
being symmetric, the broad profile requires a broadening 
mechanism other than Compton broadening, and therefore 
the relativistically-broadened interpretation is still compat- 
ible with the data. 

The inner disc radii measured from the relativistic mod- 
els range between 2 and 13 R g . We noticed some differences 
between models based on different space-time metrics. LAOR 
(for maximally rotating black holes) systematically showed 
the lowest values of Ri n , in the range 2 to 6 R g (with the 
exception of Obs. 1 where all the models measured Ri n ~11 
R g ). diskline, based on a Schwarzschild metric, showed Ri n 
between 6 and 11 R g , with one observation pegging to the 
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lowest value. Between the two models, LAOR fits the data sta- 
tistically better in most of the observations. Notice, however, 
that for rotating NS with a* < 0.3 the metric is expected not 
to de viate much from Schwarzschild metric fsee lMiller et al.l 
1 19981 . for more details), which is probably t he reason why 
LAOR under-estimated the inner disc radius. ICackett et al.l 
(|201Cf ) fitted the emission line profile on Obs. 1-3 using the 
diskline model. It is interesting to notice tha t the measure- 
ment for the inner radius and the EW from ICackett et al.l 
l |2010l ) are lower than what we reported here. On the other 
hand, line energy and emissi on index are consist ent with 
being the same. The fact that ICackett et all (|2010l ) did not 
apply corrections for pileup effects and that only used the 
XMM-Newton data to fit the Fe line, may be the reasons of 
these discrepancies with our results. 

Finally, we modelled the Fe line profile with kyrline 
that allows to set the spin parameter between (equivalent 
to a Schwarzschild metric) and 1 (maximally rotating Kerr 
black hole) . We tested three different spin parameter values, 
a, = 0, 0.27 and 1 (see Section [3.1.11 for details). With these 
models we measured values of the inner disc radius mostly 
between 5 and 13 R g (with the only exception of Obs. 5 
where for kyrline with a* =1 we found the radius to be at 
~2.5 R g )- We found consistent measurements of the inner 
radius for all three spin parameters, with only two exception 
in Obs. 2 and 5, where one of the models significantly de- 
viated from the other two (see Table [6}. It is interesting to 
mention that kyrline with a t —1 and LAOR, that are based 
on the same Kerr metric showed in general s ignificantly dif- 
feren t values for the Fe line parameters (see ISvoboda et all 
12009] , for a detail comparison between the models). Very 
similar to that, also for kyrline with a* =0 and diskline 
we noticed discrepancies in the modelling of the line profile, 
although less significant. 

None of the line parameters from the relativistic models 
show a clear trend with the source state. However, the inner 
radius of the disc appears to significantly change for two 
of the soft-state observations (Obs. 4 and 5), the two with 
the largest S z values. The line flux increased going from the 
transitional-state to the soft-state. The inner radius inferred 
from the reflection model, which is consistent with Ri n from 
kyrline with a* =0.27, did not show any clear correlation 
with the source state. 

Our findings, both from reflection and phenomenologi- 
cal models, suggest that at least for 4U 1636-53 the prop- 
erties of the Fe emission line, assuming the line profile is 
broadened by relativistic effects, are not consistent with the 
scenario in which the inner truncation radius of the accretion 
disc changes with the source state. 



4.1 Caveats of the analysis 

In all cases, the model we used to fit the six XMM- 
Newton/RXTE observations of 4U 1636-53, including rel- 
ativistically smeared ionised reflection, described the data 
well, although there are a few caveats. Firstly, we generally 
found ionisation parameter values consistent with the pres- 
ence of moderately or highly ionised Fe (log£ between 1.2 
and 3) , however the line energy values from the phenomeno- 
logical models were mostly at 6.4 keV consistent with neu - 
tral Fe. This has been also observed bv lCackett et all l|201fj ) 
between the reflection model and diskline. The combina- 




Figure 7. Transitional-state spectrum (top, Obs. 1) and two soft- 
state spectra (middle, Obs. 3 and bottom, Obs. 4) of 4U 1636-53. 
Each plot shows the simultaneously fitted XMM-Newton/RXTE 
spectra and unfolded model in the main panel, and the residuals 
in sigma units in the sub panel 

. The continuum emission components arc represented as in 
Figure [3] The relativistically smeared Fe emission line is fitted 
with the phcnomcnological model kyrline with spin parameter 
a, = 0.27 (violct/dashcd-triple dotted). 
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tion of the thermal emission from the disc and the incident 
emission are consistent with ionisation parameters we found, 
raising then doubts on the reliability of the line energy val- 
ues from the phenomenological models. However, we should 
keep in mind that with phenomenological models we over- 
simplify the reflection spectrum between 4-9 keV by only 
fitting a emission line while in reality the spectrum is a su- 
perposition of a continuum plus a line spectrum. 

Secondly, from all our fits to 4U 1636-53, we found a 
high disc inclination, which is at odds with the fact that no 
eclipses are detected in this source. From pre vious analyses 
of the first three observa t ions i n our sample, iPandel et al.l 
|200ij ) and lCackett et all (|2010T ) also reported similar incli- 
nations, always with an upper confidence li mit pegged at 
90°. T he high inclination in t he analysis of IPandel et al.l 
(|2008l ) and lCackett et all (|2010T l maybe partly due to pileup 
effects, as pileup hardens the spectrum and possib ly extends 
the blue wing of the Fe line. IPandel et al.l (|2008T ) also sug- 
gested that the broad line could be due to line blending of 
different ionisation stages of Fe. When they fit the Fe line 
with a multiple line model, they found the inclination to 
be slightly lower (> 64°), although it was still consistent 
with 90°. Nevertheless, even with a conservative pileup cor- 
rection, and fitting a model that self-consistently includes 
line blending effects and Compton br oadening (the latte r 
was not included in the model used by IPandel et al.ll2008h . 
we still found very high inclination values. Although fits to 
PN timing mode data of other NS LMXBs have revealed 
well constrained and much lower in clination angles (see e.g. 
iD'Ai et all 120091 : llaria et al.ll2009h . it is possible that the 
PN timing mode data we have used here are still affected by 
calibration issues. It would be useful to perform a spectral 
analysis similar to the one presented in this work, but ap- 
plying it to data taken in different data modes or with other 
X-ray missions. 



tic models did not consistently modelled the Fe line; this 
is especially true for the LAOR model. Values of the inner 
radius estimated form the reflection models are consistent 
with the results from kyrline with a, = 0.27. However, the 
reflection model did not help solving the issue of the too 
high inclination values derived from the Fe line profile. 
Finally, we also explored the variation of the direct contin- 
uum emission and the properties of the Fe line as a func- 
tion of the source state. According to the standard accre- 
tion disc model, as mass accretion rate increases the disc 
moves inwards. We found that the continuum emission over- 
all matched consistently with that picture, however the in- 
ner disc radius inferred from the Fe line did not show any 
significant correlation with the source state. 
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5 SUMMARY 

We analysed six XMM-Newton/RXTE observations of the 
NS LMXB 4U 1636-53; in two of the observations the source 
was in the transitional state, and in the other four the source 
was in the soft state. We used a relativistically smeared 
ionised reflection model to study the Fe emission line appar- 
ent in all six spectra, which we used to investigate the rela- 
tive contribution of the corona and the NS surface/boundary 
layer to the reflection continuum. We found that the NS sur- 
face/boundary is the main source of irradiation in the Fe line 
region in four out six observations, both in the transitional 
(also known as island state) and soft state. In the other 
two observations (one in the transitional and the other in 
the soft state), the whole reflected spectrum is due to the 
corona illuminating the accretion disc. 

We also fitted the Fe emission line with a set of phenomeno- 
logical models including symmetric as well as relativistically- 
broadened profiles. We found that, although the relativistic 
line models fitted well the data, the symmetric GAUSSIAN 
profile gave the best fit in statistical terms. The best-fitting 
Gaussian line was, however, in all cases very broad. Our re- 
sults suggested that such a broad line profile is unlikely to 
be produced by Compton broadening, and that relativistic 
broadening may still be necessary. We found that relativis- 
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Table 4. Obs. 1-3 - Phenomenological models parameters: Continuum components 
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+ 1.1 
0.9 
+ 0.004 
0.007 
+0.2 
0.1 



15.4 
0.199 
2.0 



1.85+0.01 

,+ 11 
-1.0 
n iqa+0-005 

u.iyo_ Q 0Q2 
n+0.8 

-0.6 



15. 3"* 



2.03 



1.86+0.01 
l5.5±}; a 
0.200+0.004 
2.0+0.7 



X u bC I") 1-08 (306/284) 1.18 (333/282) 1.09 (307/282) 

Total flux (10~ 9 ) 2 - 2 lo'2 2.2+0.4 



, + 0.2 
-0.1 



1.11 (312/282) 1.11 (312/282) 



9 O+0.3 
^-0.1 



9 9+0.2 
2 ' 2 -0.1 



1.12 (317/282) 
2.2+0.1 



Model comp Parameter 



Gaussian 



DlSKLINE 



Laor 



Kyrline 



Kyrline 0.27 



KY'RLINE 1 



PHABS 
DISKBB 



JV H (10 22 cm- 2 ) 
kT in (kcV) 
Afdbb 

F d (10- 11 ) 



65+ 03 
u -»5_ os 

107+ 2 ? 
— 15 
qo o + 1 1 ■ 5 

oo.o_ 15 6 



0.32+0.02 
0.67+0.02 

188 4 
66. 7 4 



„+24 



7+12.3 
-11.7 



0.31+0.03 
0.65+0.04 

127 4 
39. 6 4 



7+22 



. + 13.4 
-10.4 



0.31+0.03 
0.68+0.04 

159+ 22 
60 1 + 18 - 6 



0.31+0.03 

67+ - 02 
UD ' -0.04 
9 + 28 
-23 
.,+ 12. 9 
-14.0 



172^ 



61. 0"' 



32+ 02 
u - a2 -0.03 

67+ ' 05 
UD -0.03 

173±™ 

61 4+ 22s » 
01 -^-13.6 



kT BB (keV) 
AfBB (10- 3 ) 



1.76+0.05 
„ o + O.S 

3-»_ .5 

3 9 + 0.7 
' 5 ' 2 -0.4 



1.64+0.03 
6 l+°' 6 
5 l + °- 5 



1.75+0.05 

4 4+ ' 1 
^■^-0.2 

3.7+0.2 



1.68+0.05 
5.3+0.7 
4.5+0.6 



l.(iS~ r 



3 + 0.04 
-0.05 
5 + 06 
J ' 8 -0.5 
4.9+0.5 



+ 0.04 
-0.03 
+0.6 
-0.5 
4.9+0.5 



5. 



r 

fcT c (kcV) 
iV NTH 

Fnth (10~ 9 ) 



2.3+0.1 

6 1+ 0S 

29+ 06 
u ' 2a -0.02 
-,+0.4 
-0.3 



2.1 1 



2.1"" 



i+O.l 
-0.2 
5.4+0.2 
n+0.01 
-0.02 



0.193 



1.7 



+ 0.2 
0.5 



0.27 
2.0 



+ 0.3 
-0.2 
+ 0.06 
0.03 
+ 0.5 
0.3 



2.2+0.2 

5 6+ 1 ' 4 
°'°-0.4 

21+ 02 
U ' Z1 -0.06 
j 7 +0.3 



, ,+0.2 
^-0.3 
5 6+ ' 3 

91+ - 03 
U -^-0.05 
1 7 +0.4 
1,1 -0.8 



2.2^ 



,+0.2 
-0.1 
5 6+°' 9 
„ 21 + 0.04 
U -' !1 -0.03 
1 7+°' 4 



Xt (X 2 /V) 

Total flux (10~ 9 ) 



1.05 (298/284) 

9 »+°- 5 
2 ' 8 -0.4 



1.07 (304/282) 
9 s+°- 2 



1.06 (299/282) 

9 c + 0.5 
2 ' 8 -0.3 



1.07 (301/282) 
9 s+°- 4 



1.07 (303/282) 

9 K + 04 



1.07 (302/282) 

9 K + 05 
2 ' 8 -0.3 



Model comp Parameter 



Gaussian 



DlSKLINE 



Laor 



KY'RLINE 



KY'RLINE 0.27 



Kyrline 1 



PHABS 
DISKBB 



N a (10 22 cm' 
kT in (kcV) 
N dbb 

F d (10- 9 ) 



0.30+0.01 

79+ ' 02 
u-/ -0.08 

2+1° 
,+0.5 
-0.4 



162^ 
1.2+ 



n ,,,+0.04 
°' 36 -0.08 
0.73+0.01 

237± 3 
1.2+0.1 



0.36+0.01 
0.71+0.03 
244+51 
1.1+0.3 



0.36+0.01 
„ -3+0.03 
u ' '^-O.Ol 
24513 
1 9+O.2 
-0.1 



0.37+0.09 

„ 73+0.03 

253+38 
, + 0.3 
-0.2 



n ,,, + 0.02 
°' 36 -0.03 
73+ - 05 

244 + 



1.2 



+ 18 
23 
+ 0.4 
0.2 



kT BB (keV) 

N BB (10~ 3 ) 
F b (10~ 10 ) 



feT e (keV) 
A^NTH 

Fnth (1°~ 



1.43 



+ 0.10 
0.43 
6.9+0.2 

5.8+1.7 

, oc+0.05 
- L ' oo -0.02 

3.3+0.1 
0.15+0.01 

1.7+0.2 



1.48+0.03 
q q+0.1 
y ' -0.4 

!.3 4 ~ " 



,+0.1 
-0.3 



, + 0.03 
-0.02 
3.4+0.1 
„ ,3+0.04 
uiJ -0.02 
- + 0.5 
1 -0.2 



1.82^ 



l.B" 1 



1.33+0.07 
a 9+0.5 
8 ' 2 -2.1 
6.9+ 



+0.4 



1.82"* 



-1. 

+ 0.01 
0.02 
3.2+0.1 
+ 0.03 
0.01 



0.15 



1.6 



+ 0.3 
0.1 



1 46+ ' 06 
J "™-0.03 

9 8+ ' 1 

, Rn +0.10 
4 »0_ 03 

3.3+0.1 

0.13+0.01 
c+0.2 
-0.1 



1.5^ 



1.46 
10.4 
8.8+0 



+ 0.05 
-0.03 
+ 0.1 
0.1 



1 7 +0.2 
1 • ' — . 1 
3.3+0.1 
0.11+0.01 

1 5+ 04 
^•"-0.3 



B.7+S;g 



1.8 



a+0.2 
-0.1 
„ 3+0.9 
■ 3 - 3 -0.5 

11 + ' 04 
u " L1 -0.03 
, 3+O.6 
1 ' J -0.4 



Xt (x Z /») 
Total flux (10' 



1.09 (306/281) 
3.4+0.5 



1.19 (333/279) 

3 g + 0.5 

— 0.2 



1.12 (314/279) 

3 g + 0.5 

■ 3 '°-0.4 



1.16 (325/279) 

„ -+0.3 
o.o_ 2 



1.17 (326/279) 

„ K+0.5 
JO -0.4 



1.17 (328/279) 
3 5 +0.7 



NOTES.- A * means that the error was pegged at the hard limit of the parameter range. All uncertainties are given at 90% confidence 
level. iV dbh is defined as (Ri n / D10) 2 cos 9, with Ri n in km, Diq the distance in 10 kpc, and 8 the inclination angle of the disc. Wbb is 
Lsq/D^q, where Z/3g is the luminosity in units of 10 39 erg s — 1 . Nnth is in units of photons keV - 1 cm -2 s — 1 at 1 keV. All fluxes represent 
the unabsorbed 0.5—130 keV flux in units of erg cm -2 s — 1 . 
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Table 5. Obs. 4-6 - Phenomenological models parameters: Continuum components 





Model comp 


Parameter 


Gaussian 


DlSKLINE 


Laor 


Kyr.line 


IvYRLINE 0.27 


Kyrline 1 




PHABS 


N a (10 22 cm" 2 ) 


0.30±0.01 


0.30+0.01 


0.31+0.02 


0.30+0.03 


0.30+0.01 


0.30+0.01 




DISKBB 


kT- m (kcV) 


0.76±0.03 


0.77+0.02 


0.74+0.02 


0.78+0.02 




0.78+0.01 






N dbb 

F d (10- 10 ) 


149+ 11 

3 4+ - 8 
3 '*— 0.7 


14311 

r + 0.6 

"■°-a.5 


13611° 
2.7±0.5 


1501J 
4 +0 - 6 


153l 8 
4 1 + 13 


151+8 
4.0+0.4 




BBODY 


fer BB (keV) 


1.67+0.03 


, V o+0.02 
x - '^-O.OS 


1.72±0.04 


1.70+0.02 


1.70+0.02 


1.71+0.02 






AT BB (10- 3 ) 




9.3+1.3 


7.8±0.9 


9.6+1.6 


10.1 + 1.3 


10 1+ 13 






F b (10 _1 °) 


r o + l 1 

5 - 8 Ii.s 


6.2+0.9 


r- o + O 6 

5 - 2 Io.7 


6.4+1.1 


6.8+0.9 


^ ^+0 9 
6 ' 7 -2.1 




NTH COMP 


r 


1 9+°- 3 
0.1 


2.0+0.2 


2 l+°- 4 
" ! ' i — 0.2 


1 9+ 04 


1.8+0.3 


1 K + 0-3 






fcT (kcV) 


3.S_ 2 


4.2+0.2 


4.1+0.2 


4.0+0.2 


o n+0.3 
J ' a — 0.1 


4 0+ 1 - 7 






iV NTH 


0.10±0.02 


0.09+0.01 


14+ 003 
ul *-0.02 

6 9+ 23 


0.07+0.01 


0.06+0.01 


0.06+0.02 






Fnth (10^ 10 ) 


6 l + 1!5 


5 0+°' 9 


"111 


4 2 +13 








xl U 2 /v) 


1.08 (304/281) 


1.06 (297/279) 


1.09 (304/279) 


1.06 (295/279) 


1.06 (296/279) 


1.06 (296/279) 






lotal Hux (1U ) 


1 ^+<>- 2 
0.3 


1 r + 0.1 
— 0.3 


1.5+0.2 


1.5+0.2 


i r+0.2 
1 ' J -0.3 


1 c + 0-2 
1 '°-0.3 


Obs. 5 


Model comp 


Parameter 


Gaussian 


DlSKLINE 


Laor 


Kyrline 


Kyrline 0.27 


Kyrline 1 




PHABS 




0.28 + 0.02 


°' 30 -0.03 


n OQ +0.05 
u -^ a -0.01 


U.Z9 + U.U2 


0.29+0.01 


u.^y±u.ui 




DISKBB 


A:T in (kcV) 
A^dbb 

F d (10- 10 ) 


0.75±0.04 

71+ 19 
,1 -29 

1 5+°- 6 


0.75+0.01 

14 6lIo 
3.1+0.3 


0.73+0.02 
120+^ 
6 l +1 -° 


0.79+0.02 
146+1 1 
10.411; 4 


79+ - 01 
1541^° 
10.1 + 1.0 


0.79+g-« 
10.5±?; 2 




BBODY 


K-Ibb (keV) 


1 . y * ± u . 1 


1 £ZQ _l_ A AO 


, K _+0.04 


i .D3 + U.U3 


1.D3+0.02 


1 A _L_A AO 






AT BB (10- 3 ) 
F b (10- 10 ) 


4 9+ 06 
3 0+°- 4 


6 7+ l fi 

4 5 +1 -° 
^■°-o.i 


4.4+0.7 
3.7+0.6 


7.4+1.1 
6.2+0.9 


7.7+1.3 
6.5 + 1.1 


7 4+40 
'■*-0.4 
« 9 + 3.3 
°' 2 -0.3 




NTH COMP 


r 


2.4±0.2 


2.0+0.1 


2.1+0.1 


1 9+ 01 
x - a — 0.2 

tr+0.2 
■ 3 -° — 0.7 

0.12+0.01 

1 Q+0.2 
1,J -0.4 


1.8+0.1 


1.9+0.1 






fcT c (kcV) 
Fnth (10~ 9 ) 


4 2 +os> 
• —0.7 

29+ ' 11 

1 1+0.5 

1 ' 1 — 0.3 


3 6+ 01 
n 1S +0.04 
U1B -0.02 

1.0+0.2 


3 4+°- 8 
°- -0.1 

23+ 01 
"^-0.03 

1 Q + 02 

0.3 


3.4+0.1 

n+ 03 
uii -0.01 
i q + 0.5 
1 ' J -0.2 


3.5+0.1 
0.12+0.01 
1.3+0.2 






xl (x 2 /») 


1.12 (304/272) 


1.22 (331/270) 


1.15 (310/270) 


1.17 (317/270) 


1.18 (320/270) 


1.21 (327/270) 






lotal Hux (1U J 


n +0.5 
3.0_ Q 3 


, n +0.3 


3 n+0.2 

3.U_ Q 3 


„ „+0.2 
3 '°-0.5 


„ n + 0.5 
3 '°-0.3 


„ n +0.5 
3 '°-0.3 


Obs. 6 


Model comp 


Par tirm 't. i ■ 1 


Caussian 


DlSKLINE 


Laor 


Kyrline 


Kyrline 0.27 


Kyrline 1 




PHABS 
DISKBB 


iVpr (10 22 cm~ 2 ) 
kT in (kcV) 
^dbb 

F d (10- 10 ) 


0.32+0.05 
0.32±0.05 

o.sioi 


n ,,.,+0.03 
°' 27 -0.09 
41 + 0.04 
u -0.02 
360l 2 f 
1 5+ 1 - 1 


n o +0.05 
u ""-0.09 
n ,,+0 06 
0-33_ 05 

644+ 43B 
°* -124 

1 o +1 - 3 
-•-■"-o.e 


n qi +0.03 
u - J1 -0.06 

°- 3 4lo:°03 
611±« T 
1 1+ 1 ' 2 


n ., 1 +0.04 

u,31 -0.08 

35+ 06 
U - 3S -0.03 

581+ 422 

1 1+ 1 ' 2 


°" 29 -0.04 
0.37±0.05 
466lf/ 




BBODY 


fcT BB (kcV) 


1.79±0.13 


, 7 „+0.04 
L - ' J -0.12 


1.90+0.11 


1.89+0.08 


1.88+0.02 


1.86+0.08 






AfBB (10- 3 ) 


1.7±0.3 


, c + 0.3 


1.8+0.2 


2.2+0.1 


2.2+0.2 


2.2+0.2 






F b (10- 10 ) 


1.4±0.2 


,+0.3 


1.5+0.2 


1.9+0.2 


1.9+0.2 


1.9+0.2 




NTHCOMP 


r 


1.9±0.1 


1.9 + 0.1 


^lo.r 
16.8+$;° 


1.9+0.1 


1.9+0.1 


1.9+0.1 






kT Q (kcV) 


15 7+ 4 ' 2 
— 2.6 


14.5+1.8 


is-slt.o 


16.5+ 4 - 7 


16. 2+ 3 ; 8 








0.23±0.02 


0.19+0.01 


0.23+0.02 


0.22+0.02 


22+ - 01 
UZl! -0.03 


0.21+0.01 






Fnth (10~ 9 ) 


, 0+0.5 


2 l+°' 4 


2 4+°' 8 
^ -0.7 


2.3i: 9 : 5 7 


, o+0.4 
z,J -0.6 






xl (X 2 /") 


0.93 (259/278) 


1.06 (291/276) 


0.94 (260/276) 


0.97 (268/276) 


0.97 (268/276) 


0.98 (270/276) 






Total flux (10 -9 ) 


2 6+ - 5 


2 -6lo.5 


2 6+°- 8 


2 6+°- 5 
^'°-0.7 


2 6+°- 5 


2 6+°- 5 


NOTES. 


-A * means 


that the error was 


pegged at the hard limit of the parameter rang 


;e. All uncertainties are given 


at 90% confidence 


level. -/Vjbb is defined 


as {Ri n /Dw) 2 cos 


9, with Ri n in 


km, D10 the distance in 10 kpc. 


and 9 the inclination angle of the disc. ./Vbb is 




, where L39 is the luminosity in 


units of 10 39 erj 


5 s 1 . JVnth is in 


units of photons keV 1 cm 2 


s 1 at 1 kcV. All fluxes represent 



the unabsorbed 0.5-130 keV flux in units of erg cm 
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Table 6. Phenomenological model parameters: Fe line 







Model 


E line (keV) 


<j (keV) 


incl(°) 




3 


Norm (10 - - 1 ) 


Flux (10 -11 ) 


EW (keV) 


x'/dof 






GAUSS 


6 40+ 07 
D, ^-o.o* 


, , 7 +0.10 


_ 


- 




4 2 + os 


4 3+- 8 





C7+O.36 
°'-0.31 


306/284 






DISKLINE 


6 4o+ - 02 




9o.o±; 5 °- 






2.1 + 0.2 


2.2+0.2 





9 fi + 0. 06 
zo -0.08 


333/282 






LAOR 


6 44+° 08 
D -0.04* 




86.31-J 


10.&t°2% 


■ 

*■ -0.8 


2 5 + °- 4 


2 5 + 04 
^' s -0.3 





oc + 0.08 
ao -0.10 


307/282 


Obs. 


1 


Ky a,=0 


6 40+ 06 
D, ^-0.0* 


_ 


86.1+- 7 


10.8+. 2 ; 9 


, o + 0.4 


2.4+0.3 


9 c + 0.3 
^'°-0.5 





qo + 0.09 


312/282 






Ky a* =0.27 


6 40+ 06 
-o.o* 


- 


8S.9+- 7 


10.6+i; 9 


O , + 0.2 


2.4+0.3 


2 5 + - 4 
^•"-0.3 





00 + O.O9 
Ji> -0.07 


312/282 






Ky a*— 1 


6 40 +0.06 
-o.o* 


- 


85.1+-J 


9 - 9 ±i-.i 


,+0.8 
J± — 0.2 


2.5 + 0.3 


2.6 + 0.3 





34+0.10 
^ -0.07 


317/282 






GAUSS 


6 40+- 07 
D, ^-o.o* 


1.3±0.1 


_ 


_ 


_ 


76 + 1.1 







o3 -0.26 


298/284 






DISKLINE 


6 40+ 07 

-0.0* 


_ 


72 5+ 17 - 5 * 
' 2 '°-11.7 


io.7+*;» 


-2.7+0.2 


3.2 + 0.5 


3.3+0.5 





,,: + 0.12 
Z0 -0.07 


304/282 






LAOR 


6 40+ ' 24 
D, ^ u -o.o* 




86.S+. 3 ; 2 * 


4 0+ 50 

^-0.8 


4 0+ 14 
-0.7 


6.4+0.1 


6.5+0.3 





00 + 0.25 
Ji> -0.21 


299/282 


Obs. 


2 


Ky a,=0 


6 40+- 2 
D, ^-o.o* 






6 3 + 11 
D,J -0.3* 


i+0.5 


- + 0.9 

°- a -o.r 


r O + 0.9 
a " 5 -0.7 





42 + 0.21 
^^-0.11 


301/282 






Ky a, =0.27 


6 40+ 01 
D, ^ u -o.o* 


- 


83.9+. 2 g 


12 5+ 2 - 9 


„ + 0.7 
J ' U -0.4 


3 6+ 20 


7+2.I 
J ''-0.5 


Q 


9O+0.12 
' !o -0.09 


303/282 






Ky a* = l 


6 40+ 01 
D, ^ u -o.o* 


- 


83.0+. 21 , 


11.8lj' 2 


2 9+ 05 
2,y -0.3 


3 6+ 06 


3 7+ 06 
J ''-0.3 





,0 + 0.11 

''"-o.oo 


302/282 






GAUSS 


6 40+ * 06 

-0.0* 


1.4±0.1 


_ 


_ 


_ 


11.7+0.7 


12.0 + 0.7 





e8 -0.42 


306/281 






DISKLINE 


D, ^ u -o.o* 




90 0+°'°* 


8 4+ 07 
°- -1.5 


-2.7+0.8 


4.6+0.4 


4 9+- 4 
^ ,y -0.5 





,7+0.14 
Z '-0.06 


333/279 






LAOR 


6 40+- 05 

Di ^ U -0.0* 




87.5±=;|* 


, o + 0.2 
2 ' d -0.5 


3.5 + 0.2 


9 6+ 05 
y -°-1.0 


w.itli 





60 +0.55 
DU -0.12 


314/279 


Obs. 


3 


Ky a,=0 


6 67+ 06 

°'° -0.07 


_ 


m.2+_° % 




fi +0.5 


4.4+0.4 


' -0.5 





,q+0.15 
Z -0.05 


325/279 






Ky a* =0.27 


6.68±0.07 




86.1±°0. 7 


13.1+J- 


, + 0.4 


4.3 + 0.4 


4 6+ 04 


Q 


,0 + 0.18 
2S -0.06 


326/279 






K\' a* = l 


6.70±0.1 


- 


85.6l 2 ;l 


12. 9+ 3 ' 2 


3 4+- 5 
° -0.3 


4 s + 0-6 
^-0.3 


4 6+ 01 





,7+0.18 
Z '-0.06 


328/279 






GAUSS 


6 - 8f >_0.10 


1.0±0.2 


_ 


_ 




29 +l.l 


3 1 + 12 





,o + U.19 
2ii -0.21 


304/281 






DISKLINE 


6 43 +0.05 
Dacs -0.03* 


_ 


73.4+" 


6 + 2s 
D,U -0.0* 


-2.6 + 0.1 


2 6 + 10 
<-o -0.4 


, o+l.O 


0. 


19+0.10 


297/279 






LAOR 


6 72 +0.15 
°' ,Z -0.32* 


_ 


88.i±i;J* 


2 - 8 -o.s 


+ 0.9 
"-"-0.3 


3.5 + 0.1 


., Q + 0.1 

d ' S -0.2 





,o + 0.24 
28 -0.21 


304/279 


Obs. 


4 


Ky a,=0 


6 43+ 01 
u -^° — 0.03* 


_ 


72.3+- 2 


6 0+ 19 

DiU -0.0* 


2 4+ 02 

z - 4 -o.i 


2 4+°' 6 


2.6+0.6 





10 + 0. 08 

-"-"-o.oo 


295/279 






Ky a* =0.27 


6 43+ 07 
Dacs -0.03* 


- 


72.ll 9 ; 2 


■5 7 +2.0 


2 4+ 02 


2 4 + 06 


2 5 + - 7 
2 '°-0.4 





, 7+0.06 
-0.10 


296/279 






Ky a, = l 
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6 41+ uua 
D - 4± -o.oi* 
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_ 


89 g-89.9. 
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_ 


90.01-?* 
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6.1+0.7 


6.4+0.7 
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-0.12 
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-1.0 
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■^^-o.g 
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° -0.8* 


9 7 +0.2 
z ' '-0.1 


4 5+°- 6 


4.6+0.6 





31+0.11 


320/270 
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6 40+- 03 
D, ^-o.o* 




89.0+i; 9 * 


2.6±0.1 
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5.4+0.6 


5.5+0.2 





oc+0.15 

J0 -0.12 
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1.2±0.1 
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4 1+0-U 
^■ ± -l.l 
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4± -0.21 
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6 40 +0 02 

Di ^ U -0.0* 




90.0+- 9 * 




-2.5 + 0.1 


2 4 + 03 


, , + 0.3 
4-°-0.7 
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25+0.09 
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86.4+°;f 


6.2±1.9 


+ 13 
J ' 8 -0.5 


3.7+0.6 


3.7+0.6 
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Jy -0.12 
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6 


Ky a»=0 
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D, ^ u -o.o* 




B6.8+S-.2 


ia.a+« 


+ 0.6 
"•"-0.4 


2 6+°' 3 
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Z ' D -0.2 





,7+0.08 
z '-0.06 
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86.S+.- 7 


12.2+I; 7 


., + O.6 
^-0.4 


2.5 + 0.3 


2.6 + 0.3 
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'"'-0.07 
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"^-o.a 


2.5 + 0.3 


2.5+0.3 





,c + 0.08 
Z0 -0.07 
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NOTES.— A * means that the error was pegged at the hard limit of the parameter range. All uncertainties are given at 90% confidence level. 
j3 represents the power law dependence of the disc emissivity. The normalisation of all the models is the line flux in photons cm -2 s — 1 . Line 
flux represents the unabsorbed 0.5-130 keV flux in units of erg cm~ 2 s — 1 . Inner radius is in units of gravitational radius Rg = (GM/c 2 ). 
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APPENDIX A: PILEUP CORRECTION FOR PN 
DATA 

As reported by, e.g., iNg et all |2010l ). data can suffer from 
pileup even when the average PN count rate is below the 
critical pileup level reported by the XMM-Newton User 
handbook (800 cts s" 1 ). Triggered by that we inspected the 
dataset for pileup effects. In addition to the spectra cre- 
ated for the full 17-column wide box centred on the source 
position, we extracted spectra selecting the same events as 
before but now excising the central, 1, 3, 5 and 7 columns, 
respectively. We then used the task epatplot to determine 
whether each of these spectra have significant pileup. For 
Obs. 1 (with a count rate of ~250 cts s" 1 ), we found that 
excising only the single central column was sufficient to sig- 
nificantly reduce the pileup. For the other five (> 300 cts 
s _1 ) observations, we had to excise the 3 central columns 
to correct for pileup. Unfortunately, the task epatplot does 
not provide a quantitative measure of how pileup affects 
the spectrum, nor whether the Fe emission line suffers from 
that. To investigate the effect of pileup on the Fe emission 
line, we fitted the spectra of all datasets, including the si- 
multaneous PCA and HEXTE data, using a phenomeno- 
logical model for the Fe line (for practical reasons we will 
only discuss results about Obs. 2, that represents the soft 
state). We fitted simultaneously all the spectra previously 
used for the epatplot investigation, where we excised 0, 1, 
3, 5 and 7 columns. We left the continuum spectral shape 
parameters (see Section [3. l.ll for details on the model) free 
to vary, but we linked them across different spectra, we also 
included the simultaneous PCA and HEXTE data. By do- 
ing that we forced the continuum emission to maintain the 
same spectral shape. We only left the normalisation of all 
the continuum components free to vary in all the spectra. 
We modelled the Fe emission line with the relativistic model 
kyrline assuming a* =0.27. The line parameters were free 
in all spectra. The main assumption here is that if the data 
are not affected by pileup, the continuum spectral shape will 
not vary while excising central columns from the PSF. This, 
in principle, should also apply to the Fe line profile. In Fig- 
ure |Al] we show the data/model ratio in the energy range 
3-9 keV, after setting the normalisation of the iron line to 
0. Different colour represent different excising regions, go- 
ing from none (black points) to 7 column excised (light-blue 
points), and solid lines represent the best-fit Fe line profile. 
From Figure |A"T1 we notice that the line profile significantly 
changes with the number of excised columns. More specifi- 
cally, we find that the line energy quickly drops from 6.97 to 
6.4 keV as the number of columns excised goes from to 3. 
On the other hand excising 5 and 7 columns still gives line 
energies around 6.4 keV. The inner radius shows a similar 
behaviour, as we can notice from the line profile. Figure |AT1 
shows that pileup strongly affects the Fe line. From these 
tests, we decided for this specific observation to correct for 
pileup by excising 5 central columns of the PSF. Summaris- 
ing, based on the epatplot task and this test, we decided 
to excise 3 central columns for Obs. 1, Obs. 4-6, and we 
excised 5 columns for Obs. 2 and Obs. 3. Our extraction re- 
gions for Obs. 1-3 are more conservative in comparison with 
the regions use d in previous analysis of the same data (see 
INg et al.|[20lbT) . 

To further investigate whether pileup affected our data, 
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Figure Al. Data/model ratios in the 3—9 keV energy range, 
for different extraction regions used to generate the PN spectra 
for Obs. 2. Different colour points represent data from different 
extraction regions, while coloured-solid lines represent the best-fit 
line profiles. We simultaneously fitted six PN spectra, including 
also the simultaneous PCA and HEXTE data. We modelled the 
Fe emission line using the relativistic line model kyrline with 
a* = 0.27. The residuals are represented after setting the line 
normalisation to 0. 

we fitted the full spectra (no columns excised) and the spec- 
tra where we excised a number of the central columns, leav- 
ing all fitting parameters of the model (the same model 
used for the previous test) free to vary. We found that in 
all six observations the best-fitting parameters of some of 
the spectral components, including the iron line, changed 
significantly as we excised more and more central columns 
of the PSF, up to a point in which the changes were no 
longer significant when we continued excising columns. As 
an example, for Obs. 2, we found that the blackbody emis- 
sion from the NS surface/boundary layer and the iron line 
components changed significantly when we fit the full data 
(no columns excised) and the data where we excised the 
central 5 columns. More specifically, the temperature of the 
blackbody decreases from 1.96 ± 0.02 keV when we fitted 
the spectrum of the full PSF, to 1.68 ± 0.05 keV when we 
excised the central 5 columns. At the same time the inner 
radius of the disc inferred from the iron line changed from 
28t\\ R g to 12.5 ± 2.5 R g , the line energy from 6.97±g_ 05 
keV to 6.41Q 1 keV, and the line normalisation from 0.6±0.1 
10~ 3 photons cm -2 s _1 to 3.6±0.4 10 -3 photons cm~ 2 s~\ 
Both the decrease of the black body temperature, and the 
changes of the iron line parameters while excising the cen- 
tral columns of the point spread function are consistent with 
pileup effects. 
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